Electronic friction and liquid-flow-induced voltage in nanotubes 
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A recent exciting experiment by Ghosh et al [1] reported that the flow of an ion-containing 
liquid such as water through bundles of single-walled carbon nanotubes induces a voltage in the 
nanotubes that grows logarithmically with the flow velocity no. We propose an explanation for 
this observation. Assuming that the liquid molecules nearest the nanotube form a 2D solid-like 
monolayer pinned through the adsorbed ions to the nanotubes, the monolayer sliding will occur 
by elastic loading followed by local yield (stick-slip). The drifting adsorbed ions produce a voltage 
£J5 ' hi the nanotube through electronic friction against free electrons inside the nanotube. Thermally 

excited jumps over force-biased barriers, well-known in stick-slip, can explain the logarithmic voltage 
growth with flow velocity. We estimate the short circuit current and the internal resistance of the 
nanotube voltage generator. 

PACS 72.80.Rj, 73.63.Fg 
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£H , 1. Introduction 

In a recent, remarkable experiment [1], it was observed that the flow of an ion- rich liquid such as water through 
O bundles of single-walled carbon nanotubes induces a voltage in the nanotube sample along the direction of the flow. 
Strikingly, the voltage grows logarithmically with the liquid flow velocity over nearly six decades. The voltage 
magnitude and sign depend sensitively on the ionic conductivity and on the polar nature of the liquid, suggesting 
that ions adsorbed, or located in the immediate vicinity of the carbon nanotubes must be involved in its generation. 
This work highlights the device potential for nanotubes as sensitive flow sensors; it also presents very interesting and 
■ intriguing questions about the underlying physics. 

[ As an explanation for the observed effect, Ghosh et al [1] suggested that the dominant mechanism responsible 

f^i 1 for this highly nonlinear response to the flow could involve a direct forcing of the free electrons in the nanotubes 
by the fluctuating Coulombic field of the liquid flowing past the nanotubes, through pulsating asymmetric ratchets. 
The model assumes adsorbed ions which nonetheless do not undergo any average physical drift, and needs to invoke 
ratchets for which there seems to be no independent evidence. Earlier work [2] had suggested instead that molecular 
layers of liquid coating the nanotube physically slip along its surface, and excite a phonon wind which drags the 
free carriers in the tube. However, this model would yield an induced voltage linear in the fluid flow velocity. Back 
> in the 19th century, Helmholtz [3] had proposed an electrokinetic mechanism where a voltage induced by a liquid 
flowing past a solid substrate will appear as a streaming potential involving the ions carried by fluid flow in the 
^ ■ diffuse (Debye)layer at the interface, while the mobile charge carriers in the substrate play no role. This mechanism 
too would, for small flow velocities, predict a linear voltage increase with flow velocity. Our purpose here will be to 
propose an alternative model capable of explaining the observations - and of making prediction that could be checked 
by further work - based on flow-induced ion drift. 

The experiments where performed on nanotube bundles (the nanotubes are not aligned). The nanotubes are are 
micrometer sized in length and nanometric in diameter. Fig. 1 shows the induced voltage as a function of the 
logarithm of the water (average) flow velocity vq. The solid line is a fit to the experimental data given by 



U = U c \n\ — 



where U c = 0.26 mV and v c — 1.54 x 10 7 m/s. 
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FIG. 1. The voltage U as a function of the logarithm of the (average) fluid velocity vo- The solid line is a fit of the form 
U — U c \n(vo/v c ). Adapted from [1]. 
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FIG. 2. Carbon nanotube immersed in a flowing liquid. We assume that the liquid molecules nearest the nanotube form a 
2D solid-like monolayer, pinned to the nanotube by the adsorbed ions. As the liquid flows, the solid sliding motion along the 
nanotube is of a stick-slip nature: time intervals of elastic deformation (loading phase) are followed by rapid local yield events. 
This will result in a force on the ions which increases logarithmically with the fluid flow velocity. The drifting adsorbed ions 
will produce a voltage in the nanotube through electronic friction against free electrons inside the nanotube. In the stationary 
state there is a net flow of ions adsorbing on the nanotube upstream, and a net flow of ion desorption downstream. 



2. Ion stick-slip model 

Although the reasoning is general, let us focus here on water as the fluid, and start by assuming that ions dissolved 
in it become adsorbed on the carbon nanotube surfaces, see Fig. 2. We also postulate that the water molecules in the 
first layer around a nanotube form an (incommensurate) solid-like monolayer. Recent simulations provided evidence 
for instances where that does occur [4]. This 2D solid-like water layer will experience a pinning potential on the 
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nanotube, that is mediated through the adsorbed ions, and is in turn dragged or pulled forward by the external liquid 
water flow. The solid sliding motion along the nanotube occurs by stick-slip. Time intervals of clastic deformation 
(loading phase) are followed by rapid local yield events. During the loading phase the local stress in the 2D-solid 
near an ion will increase linear with time, until a critical stress a c is reached, when a local yield (rapid rearrangement 
process) takes place. We denote the elementary solid element where most of the rearrangement occur as the "stress 
block" . 

A linear (or nearly linear) relationship between driving force and the logarithm of the induced drift velocity has 
been observed for many systems exhibiting local stick-slip motion, e.g., charge density wave and flux line systems [5], 
as well as for atomic force microscope tips moving in a substrate pinning potential [6] . These results are collectively 
understood through a standard picture involving thermally excited jumps over force-biased barriers [9]. This picture 
will apply for all systems which can be described as an effective elastic solids pinned by defects, and driven by an 
external force. The basic picture consist of elastic loading, followed by rapid, local jumps over pinning barriers; in 
these jumps local regions (stress domains) of the elastic solid move forward. When thermal activation of such a 
stick-slip process is taken into account, the average stress at yield increases logarithmically with the velocity of the 
dragging fluid. 

Assume that the (average) velocity v of the adsorbed 2D solid-like water monolayer depends linearly on the (average) 
fluid flow velocity vq, v — avo, as expected for a Newtonian liquid. During the loading phase the strain in a stress 
block will be of order u/L, where L is the linear size of a stress block and u the displacement u — vt. The local 
shear stress in the stress block will be of order a « Eu/L, where E is the elastic modulus of the 2D solid. We define 
the loading force F = Laa = Eau = ku, where a is the thickness of the 2D solid monolayer film, and the effective 
spring constant is k = Ea. We could expect E to be similar to the elastic modulus of ice, i.e., of order 5 x 10 9 Pa 
and with a « 0.3 nm we get k w 1.5 N/m. Assuming the critical strain at yield to be of order 0.1 we get the critical 
displacement u c « 1 A and the elastic barrier e = ku^/2 w 0.05 eV, i.e., in the same range of energy as a hydrogen 
bond. 

i KftW- 



FIG. 3. Mechanical model representing the stick-slip (loading by elastic deformation followed by abrupt yield) of the 2D-solid 
in the vicinity of an ion [see Fig. 2(a)]. The spring represent the elastic deformation of the 2D solid. When the spring force 
reaches a critical value (the static friction force) the block start to slip, reverting back to the pinned state when the spring force 
approximately vanishes. 

The picture presented above is represented in Fig. 3, where the spring represents the elastic deformation of the 2D 
solid. When the spring force reaches a critical value (the static friction force) the block starts to slip and will revert 
back to the pinned state when the spring force approximately vanishes. This fast slip event corresponds to a yield 
or rearrangement process of the 2D solid in the vicinity of the ion. In [8,7] we studied the relationship between the 
drive velocity v and the (average) spring force F for the model in Fig. 3 including temperature effects. If F c denotes 
the "static" friction force (the spring force necessary in order to start the slip of the block), at temperatures T > K 
thermal excitation will allow the block to depin before the spring force has reached the critical value F c . The rate of 
thermal excitation over a barrier e can be assumed to be of the standard form 

w — iye~ f3e 

where the prefactor v w 10 10 s _1 and the effective depinning barrier e' (0 < e' < e) depends on the local shear stress in 
the contact area. By considering the statistical distribution of the microscopic shear stress, the friction versus speed 
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can be calculated [8]. Three regimes are found, depending on the magnitude of the speed. We define the parameter 
v*: 



¥ V er^ (1) 



2k(3e 

In a wide intermediate range of speeds, where 1 << v/v* << e 2/3e , we get 

4/3e \v*J w 
In the extreme low velocity limit, where v/v* « 1, we have 

F = *V« 
v 

For very high velocity instead, v/v* >> exp(2/3e), we have on the other hand 

2 V 2(3ev 

We now assume that the force F will also act on the ion which is at the same time embedded in the 2D solid 
monolayer and adsorbed on the nanotube. This force will result in a slow average drift motion of the ion along the 
water flaw direction. Since F is very small, linear response should apply. The adsorbed ion drift velocity is thus 

v a = (3) 
where the mobility \i depends on the lateral ion diffusion barriers. Combining (2) and (3) gives 

v a = vl In , (4a) 

where v c — v* /a and 

< - w < 4t > 

If we write F c = ku Cl where u c is the displacement necessary in order induce the yield, we get from (1) 



Note that (5) can be written as 



2 " e = H^I 



which shows that e is very insensitive to the exact value of a, u c and v. If we assume a = 0.1, u c w 1 A and 
v w 10 10 s _1 we get with (from Fig. 1) v c = 1.54 x 10~ 7 m/s, e w 0.2 eV which is very plausible. 



3. Flow induced current and voltage: electrical properties 

Assume that there is a 2D concentration n e of conduction electrons in the carbon nanotube. We shall use here a 
semi-classical picture for the motion of the electrons. Under the ion-originated external forces the conduction electrons 
drift with the velocity v e , so that the electric current in a nanotube will be / = 2irrn e ev e , where r is the radius of the 
nanotube. The equation of motion for the drift velocity is given by 

m e v e = -m e -v e - m a ri(v e - v a )— + eE (7) 
r n e 

where r is a Drude relaxation time, r\ an electronic friction coefficient, n a the concentration of adsorbed ions, E 
the electric field in the 2D nanotube, and v a is the drift velocity of the adsorbed ions. In deriving (7) we assumed 
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that the frictional drag from the moving adsorbates to be distributed uniformly over all the conduction electrons in 
the nanotube. The electronic friction coefficient 77 can be related to the adsorbate-induced increase in the nanotube 
resistivity. [9] Thus, when v a — (no fluid flow) we get in steady state 

eTE/m e 



l + r]T(m a /m e )(n a /n e ) 

2irrn e e 2 TE/m e 
l + rjT(m a /m e )(n a /n e ) 



so that the current 



Thus, U — RI, where the resistivity 

R = Ro (l + fp-— — =Rq(1 + k) 
V m e n e J 

where i?o is the nanotube resistivity in the absence of adsorbates: 

R = Lm e /(2nrn e e 2 T) 

We consider now the full equation (7) in two different limiting cases, namely short circuit and open circuit. In the 
first case E = so that in the stationary case (1) takes the form 

= -m e -v e - m a rj(v e - v a ) — 
t n e 

so that 

(m a /m e )(n a /n e )r]Tv a 



1 + {n a /n e ){m a /m e )riT 

and thus 

2irrn a (m a /m e )'qTev a K 
I = -!■ '- = — — 2irrn e ev a 8 

1 + K 1 + K 

The second limiting case (open circuit) correspond to v e = so that, from (7), U = LE is given by 

LmaT) n a 

U = v a — KRo2irrn e ev a (9) 

e n e 

We expect n > 1 so that for the short circuit case the current / ~ 2-Krn e ev a is nearly equal to the current expected if 
the nanotube conduction electrons drift with the same velocity as the adsorbed ions. 

We note that R given above represents the internal resistance of the nanotube bundle as an EMF generator. The 
measured resistivity po of the nanotube bundles in the normal atmosphere is po = 0.02 Sim so that the resistance 
^0 = PoL/A where the length L = 1 mm and the cross section area A — 0.4 mm 2 . Thus i? = 50 0. If the electron 
flow velocity v e is measured in m/s and if the electron concentration n e = 2 x 10 27 m~ 3 we get the electron current 
in Ampere, 

/ = An e ev e w 100 v e , 

and the induced voltage in mV, 

[/ = i? / = 5x 10 6 v e . (10) 

The lowest flow velocity for which a flow induced voltage was measured was 2 x 10~ 7 m/s and in that case the induced 
voltage was about 0.2 mV. Now, this is basically the value one would obtain if the drift velocity of the electrons in 
the nanotubes equaled the (average) liquid flow velocity and if the resistivity R for the liquid-immersed nanotube 
bundles were the same as that of the bundles in the atmosphere. Thus with » e = 2x 10~ 7 m/s Eq. (10) gives 

U = 1 mV. 
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This is only five times higher than the measured voltage indicating (if as we expect R is similar to Rq) that the electron 
drift velocity in the nanotubes, in the limiting case of the lowest possible fluid flow velocity where the induced voltage 
can be observed, is of order the average fluid flow velocity. At higher fluid flow velocity drag is less perfect, and this 
is no longer the case, since the voltage increases logarithmically with the fluid flow velocity. At the highest studied 
fluid flow velocity vo ~ 1CP 2 m/s the electron drift velocity would be roughly a factor of 10 5 lower than the average 
fluid flow velocity. 

4. Consequences of the model 

Because in our stick-slip model the logarithmic voltage increase with flow velocity is due to thermally activated 
jump of the ions, biased by the force arising from the drift motion of the solid monolayer, we can anticipate first of all 
a strong temperature dependence. This is apparent in Eq. (4) above, where the strongest temperature dependence is 
derived from the thermal Boltzmann factor which occurs in the microscopic expression for the mobility /z, that should 
thus increase exponentially with temperature. In addition there is also a linear temperature prcfactor 1//3 ~ T, but 
this is in most cases negligible compared to the Boltzmann factor. There may also be a temperature dependence 
arising from a temperature dependence of the concentration of adsorbed ions. In any case, a strong temperature 
dependence of the induced voltage is predicted by our model and should provide a first possibility of experimental 
test. 

A second aspect, which strongly distinguishes our model from fluctuating ratchets, is the flow- induced concentration 
of ions from upstream to downstream. In the free fluid, every ion is neutralized by some counter-ion of opposite charge. 
When an ion, say positive, adheres upstream to the nanotube surface, it will get to a large extent neutralized by the 
image electron charge in the nanotube. The corresponding negative counter-ion in the fluid is no longer neutralized. 
Being made redundant by the positive ion adsorption, it must go somewhere else. When on the other hand, after 
having drifted along a nanotube the positive ion leaves the nanotube surface some distance downstream, it will lose 
its neutralizing image electron, and must at that point again recuperate a counter-ion. This suggests that the flow- 
induced drag of positive adsorbed ions will be accompanied by the simultaneous flow of image electron inside the 
nanotube, and by a drift of negative counter-ions inside the fluid, that will migrate from upstream, where they are 
redundant, to downstream, where they are needed. In that case we have two negative currents, one electronic inside 
the nanotube and one ionic, due to the negative counterions, against only one ionic positive current, that of the 
adsorbed ions. The end result is a net negative current, as observed. 

This outcome differs from that of models such as that by Kumar et al, where no physical drift of ions is involved. If, 
as we suggest, ionic transport is involved, then there will not generally be a spatially uniform steady state. To be sure, 
there could be a non-uniform steady state where the ion concentration (of both signs) is slightly higher downstream 
than upstream. Depending on geometrical conditions, flow velocity etc, the ions might or might not have the time to 
diffuse backwards upstream, against the flow, so as to establish such a steady state equilibrium. That suggests the 
possibility to look in principle at that ion concentration gradient as the possible signature of an ion drift mechanism, 
although it is at the moment unclear whether such a gradient could in fact be measured. 

A third consequence of our model is the presence of characteristic stick slip noise, that should be readily observable 
by frequency analysis of the induced voltage. Straightforward as this seems, it might constitute the most direct and 
simplest test of the model proposed. 

5. Summary and conclusions 

We propose that the logarithmic voltage observed in nanotubes bundles upon immersion in a flowing liquid with 
ions dissolved, could have the same thermal barrier jump origin as that observed for frictional stick-slip processes 
between solids. Once ions adsorb on the nanotubes, they should be dragged not just by flowing individual liquid 
molecules, but by some kind of solid layer of molecules which surrounds the nanotube, and which is rigid enough 
to support the stick-slip. Each drifting adsorbed ion can in turn drag along electrons inside the nanotube, through 
the very same electronic frictional force which enables a fixed adsorbed ion to cause electrical resistance inside the 
nanotube. 

By assuming this mechanism to be at work we obtained very simple formula for the nanotube generator electrical 
properties. Expressions for the internal resistance, short circuit current and open circuit voltage obtained in this 
manner seem entirely reasonable. Moreover the latter quantities increase logarithmically with the flow velocity as is 
seen experimentally. 

Various tests can be considered for our mechanism. First, the logarithmic voltage increase should be strongly 
temperature dependent. Secondly, the physical dragging of adsorbed ions should cause the ion concentrations in the 
fluid to become non-uniform, and in particular to become higher downstream. Third, a characteristic stick-slip noise 
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spectrum should arise on top of the DC flow-induced voltage. It is hoped that these preliminary simple predictions 
will stimulate further experimental work on this intriguing effect. 
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